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FOREWORD 

The purpose of this report is to present those fundamentals of classical analytic me- 
chanics which are necessary for understanding the simulation of space vehicle motion 

pi and performance. Emphasis is placed on concepts, principles, and methods which are 

especially suitable for use in developing mathematical models of missile trajectory 

B problems for solution on digital computers. 

p For training purposes, there is presented a unified, yet detailed, discussion of coord- 

fel inate geometry (especially translations and rotations) , elementary matrix and vector 

operations, and the classical Newtonian physics governing the dynamics of particles 
y and rigid bodies. Such an integrated approach, which is required for familiarization 

of junior engineers and computer programmers with trajectory simulation, is not 

generally found in the usual textbooks on physics, calculus, and analytic geometry. 

Furthermore, this manual provides complete development of certain vital theorems 



t | which are either omitted in the textbooks or treated as obscure exercises. All topics 

Li 

are developed with mathematical rigor so that the document can serve as a source book 
ri on fundamental questions and principles employed in LMSC digital computer trajectory 

L' programs. 

n 

U This report can serve as an instructional manual for people who have studied analytic 

geometry and elementary calculus but no classical dynamics. It also can be used for 

i I ^ review study by individuals who have taken academic courses in mechanics, 

u 

tl 
ti 
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Section 1 
EULER ANGLE COORDINATE TRANSFORMATIONS 

1.1 INTRODUCTION 

Describing the full dynamics of a missile trajectory involves the use of a number of 
coordinate reference frames in which the various physical relationships are defined 
mathematically as vector components. This requires a means of specifying the relative 
orientation of one frame of reference with respect to another and a systematic method of 
performing coordinate transformations. 

Expressing relationships in matrix and vector notation not only facilitates transforma- 
tions, but preserves much physical intuition embodied in the notation. The use of co- 
ordinate transformations makes it possible and desirable to define physical quantities 
in their most convenient frames of reference. 

This section first presents the basic matrix and vector operations required to under- 
stand coordinate transformations, then describes the use of Eulerian angles and direc- 
tion cosines in accomplishing these transformations. 

1. 2 ELEMENTARY MATREK OPERATIONS 

1, 2. 1 A matrix A is a rectangular array of elements of the form 



a,. a„„ . . . a, 
11 12 In 
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n 

u 

n 



If m - n , then matrix A is called a square matrix. 



The order of a matrix is indicated by the number of rows and columns that form its 
array. For example, matrix A is of order m by n , denoted m x n , where m 
is the number of rows and n the number of columns. If m = n , the matrix is said 
to be of order m . 



1. 2. 2 A convenient way of denoting matrix A is (a.. ) 
where 



u 



i = 1, 2, . . . , m 
j = 1, *&,..., n 

and a., indicates the element in the i row and j column. 
13 



1. 2. 3 The transpose A of matrix A is 



a ll a 21 ' 



ml 



a l2 a 22 ' • ' a m2 



a, a„ . . . a 

In 2n mn 



V 

It 

c 



D 



which is obtained from A by interchanging rows and columns. 






1. 2.4 A matrix x , consisting of one row and n columns, is called an nth order 
row vector. 



t 



Y = / V Y 
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1,2.5 The transpose, x , of x is a column vector. 



x T « 



1. 2. 6 Matrix Sum and Difference. The sum (or difference) of two m x n matrices 
A and B is defined as follows : 



A ± B = 



11 



12 In 



a 21 a 22 



l 2n 



8. - £L « 

ml m2 



mn 



b b 

11 12 



b b 

21 22 



ml m2 



In 



2n 



. b 



mn 



a il * b ll a 12 ± b 12 



a 21 * b 21 a 22 * b 22 



• • a m * b m 



. . a„ ± b„ 
2n 2n 



a , ± b , a ± b „ . . . a ±b 
ml ml m2 m2 mn mn 



EXAMPLE: 
If 



r 
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then 



A + B 



2 3 

1 -4 

-3 12 



12 

4 0-1 

3-2 3 



3 2 3 

5 -4 -1 

0-15 



1. 2. 7 It is evident from an inspection of paragraph 1. 2. 6 that addition of matrices is 
commutative and associative; that is, A+B = B+A and A + (B + C) = (A + B) + C 



1. 2. 8 Matrix Products. The product of two matrices A and B is defined only if the 
number of columns in A is equal to the number of rows in B . If A is an m x n 
matrix and B is an n x r matrix, then their product AB is an m x r matrix defined 
as follows: 



AB = 



a ll a 12 



a 21 a 22 



ml in2 



a 



In 



a 



2n 



a_ 



mn 



b ll b 12 



b b 

21 22 



b . b „ 
nl n2 



a^, b n , +a in b„ 1 + . . . + a, b , 
11 11 12 21 In nl 



a 21 b ll + a 22 b 21 + • • * + a 2n b nl 



a ,b,- + a chni + • 
ml 11 m2 21 



+ a b - 
mn nl 



Ir 



2r 



. b 



nr 



a ll b 12 + a 12 b 22 + ' * • + a ln b n2 



a 21 b 12 +a 22 b 22 + • • • +a 



ml 12 m2 22 



2n b n2 



. +a b 

mn n2 



u 



T i 



D 



t-v 
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EXAMPLE: 
If 

A = 



2 3 

1-4 

-3 12 



then 



AB 



and 



AC = 



2 3 

1-4 

-3 1 2 



B = 



12 

4 0-1 

3-2 3 



12 

4 0-1 

3-2 3 



, C = 



1 
4 
3 



11 -2 9 

-15 2 4 

7 -10 5 



2 3 

1-4 

-3 12 





1 




11 




4 


= 


-15 




3 




7 



1. 2. 9 On inspecting paragraph 1.2.8, we see that multiplication of matrices is associ- 
ative but not, in general, commutative; that is, (AB)C = A(BC) but AB ^ BA , 
in general. 

1. 2. 10 Further inspection reveals (paragraphs 1. 2. 6 and 1. 2. 8) that matrices are 
left and right distributive with respect to addition; that is , 
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f * 



U 



1.2.11 If is the matrix (b..) such that b,. = for all i and j , then is 
called a null or zero matrix and A + = A . 






1. 2. 12 The expression I = (c. ) is called an identity (or unit) matrix when it is a 
square matrix whose elements are all zeros except for those on the principal diagonal 
(top left to lower right) which are all ones : 



1 . . . 
1 . . . 
1 ... 



(c.) = 



. 



n 

n 
11 



It can be shown that AI = A and IB - B , where the indicated products are defined. 
If A is a square matrix, then AI = IA = A . 

-1 -1 

1. 2. 13 The inverse of matrix A is denoted A and has the property that both A A 

_1 

and AA are unit matrices. 

T 
1. 2. 14 A matrix A is orthogonal when AA = I , an identity matrix. If A is orthog- 

T 
onal, then its transpose is equal to its inverse; that is, for A orthogonal, AA = I , 

A _1 (AA T ) = A _1 (I) , (A _1 A) A T = A -1 , IA T = A -1 , and hence A T = A -1 . 



fl 





1. 2. 15 Linear Transformations. If A is an n x n matrix and x is an n order 



column vector, then A transforms x into y as follows: 
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then the linear transformation of x by A is 



Ax = 



l ll 



a 12 a 13 



a rt „ a„„ a„„ 
21 22 23 



a 31 a 32 a 33 





— - 






x l 






X 2, 


= 




x 3 






_ • - 





a ll X l + a 12 X 2 + a l3 X 3 



a 21 X l + a 22 X 2 + a 23 X 3 



a 31 X l + a 32 X 2 + a 33 X 3 



If A is orthogonal, then y can be transformed back into x as follows: 



- -1 - -1 

Since Ax = y , then A Ax = A y , and 



x = A y = A y = 



a ll a 21 a 



a 12 a 22 



a 13 a 23 a 33 



— 




— — 




31". 




y l 




32 




y 2 


= 


33 




y 3 





a il y l + a 2l y 2 + a 31 y 3 



^#1 + a 22 y 2 + a 32 y 3 



a 13 y l + a 23 y 2 + a 33 y 3 



1.3 ELEMENTARY VECTOR OPERATIONS 

1. 3. 1 A scalar is a quantity that is characterized by magnitude only or by magnitude 
and an algebraic sign; e.g. , mass, time, distance. 

*• 3. 2 A vector is a quantity which has direction as well as magnitude; e. g. , force, 
velocity, acceleration. It is represented graphically by an arrow whose length indicates 
magnitude (or modulus) and whose terminus indicates direction. 



It follows from this definition that a vector quantity is preserved when displaced par- 



all £*1 +r\ •t + ^.-.l^ 
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Vector quantities are denoted by a letter with an overline; e.g. , a . The magnitude 
only of a vector quantity is denoted by the letter without the overline, or by indicating 
the absolute value of the vector; i. e. , a or | a| . In matrix notation, a vector is 
usually columnar (paragraph 1. 2. 5) with each element denoting the value of the respec- 
tive components of the vector in a coordinate space. 

A vector a in a right-handed system of coordinate axes OXYZ has its magnitude and 

direction completely determined by 
the three component vectors a , 
a , a in the respective directions 
of the coordinate system (Fig. 1). 

In Fig. 1, it is apparent that the 
modulus of a is 




-/ 



2,2,2 

a + a + a 
x y z 



X' 



Fig. 1 



If a , 8 , and y denote the angles between a and the three component vectors a , 
a , and a , respectively, then 













n 



e 

D 
D 



a x a v a z 

cos a - — , cos 8 = — 5L , and cos y - ~— 
a r a 'a 



These cosines are known ss the direction cosines of the vector a with respect to the 
the coordinate system OXYZ. 



p 

•1 



n 
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i- 3 - 3 Addition of two vectors a and b is shown in Fig. 2. Construct the origin of 

g one of the vectors (e. g. , b ) so that it is 

coincident with the terminus of the other. 
The resultant vector c , determined by 
constructing a vector from the origin of a 
to the terminus of b , is the sum of a and 
Fig. 2 b . 




lis! 



1 



&3 



Numerical values can be obtained by employing trigonometric identities. In this man- 
ner, any number of vectors can be summed. 



If b (Fig. 3) is the vector which when added 
to a will give c , then a = c - b is the 
difference between c and b . 




Fig. 3 



1.3.4 A vector v multiplied by a scalar s results in the vector sv , which is s 
times greater in magnitude than v and in the same direction if s is positive, but 
opposite in direction if s is negative. 



1. 3. 5 The vectors i , j , and k (Fig. 4) are unit vectors in the direction of x , y , 
and z , respectively, of OXYZ. They are defined as having the property 
Z 



6— 



K. i i 




I i| = 1 , |j I = 1 , |k| = l 

Thus, if the projections (Fig. 4) of the vector 
a on the coordinate axes are 
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a. - ia , a = ia , a = ka 
x x y J y z z 



t: 



Moreover, since a = & + a, + a , a can also be given by a = la + ja + ka . 

x y z a j x J y z 



1-3.6 Scalar (Dot) Product of Two Vector s. The scalar or dot pro duct of two vectors 
is equal to the product of their magnitudes, multiplied by the cosine of their included 

angle (see Fig. 5). 



Symbolically , 

a • b = J a. | |b J cos y 




Fig. 5 



u 

n 



n 



It follows that scalar products between unit vectors give the special relationships 



i • i = j • 3 = k • k = 1 and i- j = i*k=j*k = 







Consequently, for 



a = i a + j a + ^ anc j b = ib + jb '+ kb 
x J y z x J y z 



a ■ b = ( ia„ + 3a, + ka ) • ( ib + 3b + kb ) 
x J y z ' % x J y z' 



u 



B 



n 



= (ab ) i • i + ( a b ) I • j + ( a b )]• k + (ab ) j • i 
v x x' v x y ' J v x z ; % y x' J 

+ (a y b y ) j • ] + (a y b z ) j • k + (a z b x )k • i + (a z b y )k • j 



t • 



+ (a b )k • k 
x z z' 



The scalar square of a vector gives 
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a ■* a = 



a| cos (0) 



= a 



a x b 



1. 3. 7 Vector (Cross) Product of Two Vectors. The vector or cross product a x b 
of two vectors a and b (Fig. 6) is equal to the vector whose magnitude is the product 

of their magnitudes multiplied by the sine of 
their included angle. The cross product is 
perpendicular to the plane of the two vectors 
such that a rotation of vector a about the 
cross product into vector b would effect a 
right-hand screw to advance in the direction 
of the cross product. Symbolically, 




PLANE OF 
a AND b 
©> a ^"' 



Fig. 6 



a x b | = | a j |b J sin 8 



Since the rotation that carries b into a is opposite to that which carries a into b , 
then 



axb = - b x a 



If e - or 180° , then 



la x bl = I al I bl sin 8=0 



Also, 



- - - 2 

la x al = |a! sin 8=0 



I 1 
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Cross products between unit vectors give 



* 1 

kJ. 



i x i 


= j x j = k x k 


i x j 


= k = - j x I 


I x k 


= I = -k x j 


k x i 


= j = -I x ic 



Consequently, for 



a = i& + ya. + ka and b = ib + ib + kb 
x y z x J y z 

a x b = ( ia + ja + ka ) x ( ib + jb + kb ) 
x y z .x J y z' 

= (a x b x ) ixi + (a x b y )Ixj + (a^Ixk + (a y b x ) jxl + (a y b y )jxj 

+ (ab ) jxk + (a b )kx i + (a b )kx j +(ab )kxk 



12 



•U** 



y z' 



Z X 



z y 



z z- 



n 






= (a x b y )k - (a x b z ) j - (a y b x )k + (a y b z ) i + (a^) j - (a z b y ) i 



f'| 



= (a y b z - ^y)" 1 + < a z b x - a x b z)3 + (a x b y - a y b x )ic 
which can be written In determinant form as 



a x E 



a a 

x- y 



b b b 

x y z 



n 







F7 



f! 



-~t 



3-80-61-2 



1 



3 



1-3.8 Radius Vector. Consider a point (x , y , z ) in the right-hand system OXYZ (Fig. 

7). The vector determined by the origin O 
and the point ( x , y , z ) , where ( x , y , z ) is 
the terminus , is called the radius vector r 
of point ( x , y , z ) , and 



r = ix + jy + kz 

If r indicates the motion of a point with 
respect to time, then, in function notation, 

r = r(t) 

and x=x(t),y=y(t),z=z(t) 




Fig. 7 



1-3.9 Vector Derivatives. Consider a vector a which changes in magnitude and direc 
tion correspondingly with a change in a scalar variable t ; that is , 

a = f(t) 



If t increases infinitesimally by At , then a changes infinites imally by Aa (Fig. 8). 

Symbolically, 

a + Aa = f(t + At) . 

Hence, the vector differential of a is 




Fig. 8 



Aa = f(t + At) - f(t) 



3-80-61-2 



n 






da = lim Aa _ lim f(t + At) - f(t) 
dt At~0 At At-0 At 



which is the vector derivative of a with respect to t . 

The second derivative of a is obtained in the same way, as in 

A 



2- 2- 

d a _ lim A a _ lim 

-.2 At-0 A ,2 At-0 
at At 



f(t + At) - f(t) 
At 



At 



Higher order derivatives of a can be similarly obtained. If k is a scalar variable 
which is also a function of t , that is , 

k = F(t) , 

then the differential of the product u of the scalar function k and the vector function 
a = f(t) is 

A(ka) = Au = F<t + At) f(t + At) - F(t) f(t) 

= F(t + At) f(t + At) - F(t + At) f (t) + F(t + At) f(t) - F(t) f (t) 



n 

u 



ft 

u 

V 

L 



t: 



D 



and the vector derivative of u with respect to t is 



du 

dt 



lim Au 
At - At 



IT n F(t + At) 

At — x ' 



f(t + At) - f(t) 
At 



lim 



F(t + At) - F( t) 



At 



f(t) 
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Consider the two vectors a and b which are functions of a scalar t. The scalar 
product of a and b is 



a-b=ab. + ab + a b , 

xx y y z z 



and the vector derivative of a • b is 



d(a • b) _ x . , x 



dt 



db da db da db z da^ 



dt xdt y dt y dt z dt z dt 



db , da r 

= a • aT + ~ dJ ■ b 



The vector product of a and b is 



a x b = (a b - a b J!" + (a^ - a^)] + (a x b y - a y b x )k 



and the vector derivative of a x b , where i , j , and k are constant, is 



d(a x b) _ _d_ 
dt dt 



<\ b * ' \ b j )l + (& A " a xV ] + < a x b y " a y V\ 



= < a y b Z - ^V* + (a y b z - ; Z b y ri + < a A - a xV 3 

+ < a z b x - a xV 3 + < a x b y - a y b x> K + < a x\ " a yV R 



db , da „ r 
a x dF + dt x b 
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1. 4 EULER ANGLES AND DIRECTION COSINES 



i ■ 



1.4. 1 Consider a point P relative to the right-handed coordinate system OXYZ 

(Fig. 9). The angles a , fi , and y 
between the ray OP and the X , Y , and 
Z axes, respectively, are direction ang les 

of OP . The cosines of these angles 

I 

cos a = r , cos//5fj =•*■ and cos y = -' 

s ■ \_y s> ? 




/5 



'4 



where 



"4- 



x + y + z 



Fig. 9 



F 



i 



w 

I.- 

n 

o 



are called the direction cosines of the ray OP. 



Now consider two right-handed coordinate systems OXYZ and OX'Y'Z' which have 

a common origin (Fig. 10). If the angles 

Z 

between respective coordinate axes are 
given by 






OX 


OY 


OX' 


a l 


h 


OY' 


a 2 


^2 


OZ' 


a 3 


^3 



To 



then the direction cosines of the OX' , OY' , 



n 



r~\-%r r\-ir 
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i 

' J 

i 



X 

C , = cos a. = - 
x'x 1 s 



\ ' C x'y = C0S h = ^ • C x<z = cos >! = ^ 



C y'x = C0S a 2 = if> C y'y = cos ^2 - if' C y'z = C0S ?2 " if 
C Z 'x = COS «8 = ^ • C Z <y = C0S ^3 B if V C z.z = C0S ?3 = if 



where 



/~2~~ 2~ 2 /~2 2 i" f 

! i = \/ x i + yi +z i' s 2 =n/ x 2 + y 2 +z 2> s s B \/ : 



x 9 + y ; + z ; , s Q - Jx| +y| + z| 



and the points (x % , y x , Zj) , (x 2 , y 2 , z g ) , (x g , y g , z g ) are on the OX' , OY' , 

and OZ' axes, respectively. 

Using vector notation, the point (x^ y 1 , z ) is the terminus of the vector 

i x = x x + 7j + Zl 
where 



Xj = i x 1 , y^ = j y 1 ', and Zl = kz. 



Since s. lies along the OX' axis, the unit vector P along this axis is given by 
^ _ s. ( x 1 +• y 1 + zL 

. | it = _i = _i i__i 

S l S l 
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The above process, then, produces the three unit vectors 

1 i' = I cos a + j cos p - + k cos y 1 



f ! 



j 1 = I cos a- + j cos /3 + k cos y„ 



'fl k' = i cos a r + j cos /? + k cos y„ 

P o o o 



which give the following table of direction cosines : 



i ■ 



V 



c , c 

x'x 



x'y x'z 



y'x y'y y'z 



z'x z'x z'z 



G 



The nine direction cosines can be expressed as the matrix 



A = 



C C C 

xx xy xz 



c c c 

yx yy yz 

c c c 

zx zy zz 



where the prime superscripts have been omitted for brevity. The transpose of matrix 
A is then 









V 



Consider the matrix product 
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T 
AA = 



C C 

xx xy 



yx 



zx 



yy 



zy 



xz 



yz 



zz 



XX 



xy 



xz 



yx 



yy 



c 



yz 



zx 



zy 



(c 2 + c 2 +c 2 ) (c c +c c +c c Vc c +c c +c c \ 

^ xx xy xzj \ xx yx xy yy xz yz/\ xx zx xy zy xz zz! 

) lc 2 +C 2 +C 2 ) (C C +C C +C C \ 

/ \ yx yy yz/ \ yx zx yy zy yz zz| 

) (c 2 + c 2 + c 2 ) 

"z/ \ zx zy zz/ 



CC +CC +CC 

yx xx yy xy yz xz 



(c „ +cc +c c Vc C +CC +CC 

\ zzCxx zy xy zz xz/\ zx yx zy yy zz y: 



Each element of the above matrix can be shown to be equivalent to some scalar product 
of two of the unit vectors i' , j' , k' . For example, 



it • il = 



i' = (iC +jC + k C ) • (iC +JC +kC ) 



XX 



xy 



XZ 



XX 



xy 



xz 



= c 2 + c 2 + c 2 

xx xy xz 



By substitution, there results the matrix of scalar products 



i< . i> 



i' • j' i' • k' 
j' • j' j' • k» 



k' ■ i' k' . j ! k' • k' 
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n 



Hence, the product 



T 
AA = 



10 
10 
1 



which is the identity matrix I (Section 1. 2. 12). 

Therefore, the matrix A of direction cosines is orthogonal (Section 1. 2. 14), and its 
transpose is equal to its inverse. 

Suppose now that a vector p is given by 

I p = i'x' + j'y' + k'z' 



n 
Li 







Then, substituting the preceding relationships for I' , j' , and k' we obtain 



p = ix'C + jx'C + kx'C + iy'C + jy'C + ky'C + iz'C 

xx J xy xz J yx JJ yy J yz zx 

+ jz'C + kz'C 
J zy zz 

again omitting the superscripts from the C's for brevity. Regrouping in terms of i , 
j , and k , we obtain 









p = i(x'C + y'C + z'C ) + j(x'C + y'C + z'C ) 

v xx J yx zx ' J v xy J yy zy ' 



+ k(x'C + y'C + z'C ) 

x yz J yz zz ' 



IL.J 



y 



bh 



,J 



, J 



J 
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The quantities in parentheses are thus precisely the components of the vector p ex- 
pressed in the OXYZ coordinate system. If we call these coordinates x , y , and z , 
then 



X = X'C + y'C + z'C 

xx J yx zx 



y = x'C + y'C + z'C 



xy 



yy 



zy 



z = x'C + y'C + z'C 
yz J yz zz 



These relationships can be written in vector -matrix notation as 



x 



C C 

xx yx 



zx 



C C C 

xy yy zy 



C C C 

xz yz zz 



y' 



The matrix of direction cosines is said to be the matrix which transforms the OX'Y'Z' 
coordinate system into the OXYZ coordinate system. 

Since this matrix is orthogonal , we can write 



y f 



L" J 



c C c 

xx yx zx 

c c c 

xy yy zy 

c „ c c 

xz yz zz 



-X 


X 






y 


S= 




z 





c c c 

xx xy xz 

c c c 

yx yy yz 

c c c 

zx zy zz 



y 



C x+Cv+C z 
xx xy* 7 xz 



C x+Cv+C z 

yx yy' yz 



C x+Cv+C z 

zx zy zz 
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n 

fcssi/ 



Another method of relating the two coordinate frames OXYZ and OX'Y'Z' involves the 
use of Euler (Eulerian) angles. These angles are obtained by determining the succes- 
sive angles of rotation through which one frame must be rotated in order to make it 
coincident with the other frame. For example, OXYZ can be brought into coincidence 

with OX'Y'Z' by: 

NOTE 

A different perspective is used 
■ in Figs. llthroughl3frcm that 
used in Fig. 10. 

(1) A rotation through an angle <p shout OZ 
to make OY coincident with OV, which is 
the line of intersection of the planes 
OXY and OX'Y' (Kg. 11); 

(2) A rotation through an angle 6 about OV, 
the new position of OY, to make OZ co- 
incident with OZ' (Fig. 12); and 




PLANE OF / 
ROTATION-^ 



x 



V 



OX'Y' 

PLANE 



Fig. 11 




PLANE OF 
ROTATION 

X' — 

vj-" 




(3) A rotation through an angle ip about OZ' 
(the new position of OZ) to make OV(the 
new position of OY) coincident with OY' 
(Fig. 13). 

During the three rotations , OX went in turn 
to OW, OW , then OX' . 




n 

p 

U 



n 

iJ 



u 
n 

LJ 



r J 



vJ 
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The angles <f> , 6 , and ^ are the Euler angles. The relative orientation of two 
ordinate systems in terms of Euler angles is developed in the next section. 



co- 



rt 



'] 



1A - 2 ^olii^PiMA^an^Yaw^ Consider -a single rotation through angle <f> of axis OY' 

from axis OY (Fig. 14). The direction cosines 
Z' Z relating the two frames OXYZ and OX'Y'Z' 

would then be 



P(x,y , z) 



OX' 
OY' 
OZ' 




Fig. 14 
OX OY OZ 



C C C 
xx xy xz 

C C C 

yx yy ^y Z 

c c c 

zx zy zz_ 



cos <p 



cos (|-0) 



,7T 



COS (p+0) COS 



cos sin ^ 



-sin <p cos 



>^ 






The point P can be referenced to the new frame through the following relations of 
coordinates x' , y' , and z» : 



x 1 = xC + yC wr + zC - x 
xx J xy xz 



y i = xC ^ + yC ^ + zC ^ = y cos + z sin 
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10 

cos <p sin <p 
-sin <p cos <p 



or 



1 

cos (j> - sin <j> 
sin <p cos <f> 



x' 

y' 






n 

R 



The matrix 



10 

cos <p sin (p 

-sin <p cos 



is commonly known as the roll matrix; that is, a rotation of the YZ plane about the X 
axis is known as a roll. 

1.4.3 Generating the X'Y'Z' frame by a rotation of the XZ plane through an angle 

6 about the Y axis (Fig. 15) results in the 
direction cosines 



X 




y Fig. 15 







OX 


OY 


oz 


OX ! 




XX 


C 

xy 


c 

XZ 


OY' 




C 

yx 


C 

yy 


C 

yz 



cos 6 




IT 



COS ( -k + e ) 



cos -sin 6 



Pi 

L 



L ■ 



D 

D 
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which is commonly known as the pitch matrix, and also results in the referencing of 
point P to the new frame through the following relations of coordinates x' , y' , 
and z 1 : 



x' = xC + yC + zC = x cos 6 - z sin 
xx J xy xz 



y t = x c + yC + zC = y 
J yx J yy yz J 



»-+' 



z ! = xC + yC + zC = x sin 6 + z cos 

zx J zy zz 



These equations may be expressed as 



cos 0- sin 

10 

sin 6 cos 6 



or 



cos 6 sin 

10 

-sin 6 cos 6 



x' 

y' 

z ! 



1.4.4 Generating the X'Y'Z' frame by rotation of the XY plane through an angle $ 
Y about the Z axis (Fig. 16) results in the 

direction cosines 




Fig. 16 



I OX OY OZ 
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which is commonly known as the yaw matrix, and also results in the referencing of 
point P to the new frame through the following relations of coordinates x» , y' , 
and z' : 



n 



X ' = xC xx + yC xy + zC xz = x cos ^ + y sin ^ 



y' = xC + yC + zC 

yx * yy yz 

z' = xC + yC + zC 

zx J zy zz 

These equations may be expressed as 



-x sin ip + y cos !p 



x' 



y' 



cos if> sin # 



sin ip cos # 



0" 




x 







y 


1_ 




z 



or 



cos # -sin ip 



sin ^i cos ip 



x' 



y' 



n 
r.J 

r i 
11 



1. 4. 5 Suppose we wish to accomplish a roll, followed by a pitch, followed by a yaw. If 
A is the roll matrix, B the pitch matrix, and C the yaw matrix, then the process 
would be 



Roll 



"x 1 " 




x 


y' 


= A 


y 


_z' 




_ z _ 



Pitch 



"x M " 




x' 


y" 


= B 


y' 


z"_ 




_z'_ 



Yaw 



"x ,M " 




"x"" 


yllt 


= C 


y" 


Z nt 




z" 



n 

U 

y 



i; 



D 
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Bv substitution, we obtain 



1 



"x 1 "" 


y m 


_z'"_ 



= CBA 



y 



Thus, it would be possible to first perform the matrix multiplication CBA (which is 
not, in general, commutative), thereby resolving the roll, pitch, and yaw to the single 
matrix M , as in 



"x m 




X 


y m 


= M 


y 


_Z IM 




z 



= CBA 



x 



1. 4. 6 A roll about the X axis, followed by a pitch about the resulting Y' axis, followed 
by a yaw about the resulting Z' axis may be accomplished as follows : 









Yaw 






to 


Pitch 






~x r 




COS Ip'' 


sin tp ; 


o" 




COS f> 





-sin f~ 




y' 


= 


-sin ip 


COS tp 










1 







. z '_ 










1_ 




sin 9 





cos & 





Roll 
"l ■ 0. 

cos <p sin <p 

>' -sin cos (i> 



cos 4> 



sin ip 



cos 



sin 0'sin 6 



-sin cos (b 
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* > 



t » 



cos ip cos 
-sin 4 ! cos 
sin 



cos ip sin sin $ + sin ^ cos $ 
-sin i}> sin sin <p + cos ^ cos <p 
-cos sin <p- 



-cos ip sin cos <b + sin ^ sin cp 
sin ^ sin cos <fe + cos # sin <p 
cos cos <p 



Ij 



f! 
tj 



p 



E) 



In what follows, we show that the nine elements of the transformation matrix in terms 
of the Euler angles 4> , 4' > and are precisely the direction cosines relating the 

OXYZ and OX'Y'Z' frames. 

For simplicity, let the above transformation matrix be represented by the matrix 



D - 



D 



xx 



D 



xy 



D 



xz 



D 



yx 



D 



yy 



D 



yz 



D 



zx 



D 



zy 



D 



zz 



Let us define three vectors p 1 , p ? , and p~ with the following components in the 
OXYZ coordinate system: 



p. = ID + jD + kD 

*i XX J XV xz 



n 

u 

n 






P 

n 
IJ 

H 

Li 

P 
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We now find the components of p in the OX'Y'Z' coordinate system by applying the 
transformation D to p 



I 



Pi = 



D 



xx 



D 



xy 



D 



y* 



D 



yy 



D 



zx 



D 



zy 



D 



D 



D 



) 
xz 




XX 


yz 




D 

*y 


zz 




D 

xz 



= (D 2 + D 2 + D 2 ) I' + (D D + D D + D D ) 1' 

v xx xy xz ; v yx xx yy xy yz xz' J 

+ (DD+DD-+DD)k» 

zx xx zy xy zz xz' 

By actual substitution of the trigonometric identities for the elements of matrix D 
(which is left as an exercise for the reader), we obtain 



P 1 - V + Oj 1 + Ok' = i' 

Thus, p 1 is the unit vector along the X 1 axis of the OX'Y'Z' coordinate system. 
Similarly, we can show that p = ]' and p = k' . 



The angles between the OX'Y'Z' and QXYZ axes were listed previously and are repeated 
here for convenience. 





i 


I 


k 


i' 


a l 


h 


y l 


I' 


a 2 


h 


y 2 
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Now 



But 



1' • 1 = COS Q.' = C 

1 XX 



IP 

I 



v . 



i • i - P 2 • i - ( iD^ + jD + kD ) • ( i + 0j + Ok) = D. 



xy xz- 



XX 



D 



Therefore , 



n 



D = C 

XX XX 



Similarly , 



n 



but 



i' • j = cos £_ = C 



xy 



11 ' ] " Pi ' ] = (lD xx + ~i% + &W • (Oi + J + Ok) = D xy 



and a therefore, 







M 



D = C 

xy xy 



By similar means, we can show that 



D D D 

xx xy xz 

D D D 

yx yy yz 



c c c 

xx xy xz 



= ic c c 

yx yy yz 






n 






1 
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Thus, the elements of matrix D are the direction cosines relating the OX'Y'Z' and 
OXYZ coordinate systems , and we have show how to compute the direction cosines 
from the Euler angles <p , 6 , and i\> . 



We can now simplify our notation in expressing the roll, pitch, yaw transformation to 



x' 



y' 



c c c 

xx xy xz 



c c 

yx yy 



yz 



c c c 

zx zy zz 





X 






y 


s= 


. 


z 





C x + C y + C z 

xx xy J xz 



C x + C v + C z 

. y x yr yz 

C x + C y+C z 
zx zy J zz 



Since matrix multiplication is not, in general, commutative, it is essential that the 
transformations be accomplished in the order implied. The example above exhibits 
first a roll (by showing it as the first multiplier of the column vector), then a pitch, 
and finally a yaw. The number and order of rotations is determined by the nature of 
the physical problem. 



